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Phase structure of an Abelian two-Higgs model 
and high temperature superconductors 
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We study the phase structure of a three dimensional Abelian Higgs model with singly- and doubly- 
charged scalar fields coupled to a compact Abelian gauge field. The model is pretending to describe 
systems of strongly correlated electrons such as high-T c superconductivity in overdoped regime and 
exotic phases supporting excitations with fractionalized quantum numbers. We identify the Fermi 
liquid, the spin gap, the superconductor and the strange metallic phases in which densities and 
properties of holon and spinon vortices and monopoles are explored. The phase diagram in the 37? 
coupling space is predicted. We show that at sufficiently strong gauge coupling the spinon-pair and 
holon condensation transitions merge together and become, unexpectedly, first order. 
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The physics of high-T c superconductivity Q, Q is not 
yet completely understood. However, certain outstand- 
ing features are generally recognized. At normal temper- 
atures, all known high-T c superconductors are ceramic 
materials characterized by a poor conductivity. At low 
temperatures the clean material is a Mott insulator. Dop- 
ing it with impurities, the insulator may become a su- 
perconductor at low enough temperatures. The physics 
of high-T c superconductivity is essentially a 2D phe- 
nomenon since charge carriers - which are electrons or 
holes provided by the dopants - are confined to the CuC>2 
planes typical for all high-T c superconductors. 

A popular approach to high-Tc superconductivity is 
provided by the t-J Hamiltonian |3j describing hopping 
holes and localized spins in a plane: 



H tJ = 
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Pi 
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The first term describes electrons moving without chang- 
ing spin a. Double occupancy is explicitly forbidden by 



the presence of the projectors Py 
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0(1- 



The second term specifies the anti-ferromagnetic Heisen- 
berg coupling between spins located at the copper sites. 
Here Si = (1/2) J^aa' 4a^a'C la ' is the spin operator, c| CT , 
Cia are the hole (or electron) creation and annihilation 
operators, n^a — c\ a Ci a denotes the occupation number, 
and rii = riirf + n,-^. 

The t-J model (JTJ is often treated in the slave-boson 
technique |4J, which splits the spin and charge degrees 
of freedom of the electrons. The electron creation op- 
erators are decomposed as 0, Q c| ff = f^h, where 
f ia is a spin-particle ("spinon") operator and &j is a 
charge-particle ("holon") operator, fn order to forbid 
double occupancy of sites one imposes the constraint 
/1/if + f^ifn + b\bi = 1 on the physical states of the 
system. The constraint drastically simplifies the treat- 
ment of the t-J Hamiltonian JJJ. 



In addition to the ordinary electromagnetic (external) 
gauge symmetry, the spin-charge separation naturally in- 
troduces an (internal) compact U(X) gauge freedom, 



U(l) 



hit 



fi< 



' fit 



(2) 



which plays an essential role [2j in understanding the 
physics of strongly correlated electrons. Besides con- 
densed matter physics, the spin-charge separation idea 
is also applied to strongly interacting gluons in QCD 0]. 

The emerging effective theory of superconductivity can 
further be simplified and reformulated in terms of lattice 
gauge models 0, H S E3- Thus, the t-J model JTJ is 
related to a compact Abelian gauge model with internal 
symmetry which couples holons and spinons. As in 
usual BCS superconductivity, under appropriate condi- 
tions the spinons couple and form bosonic quasiparticles. 
In a mean field theory one can define fields which behave 
under the gauge transformations @ as: 



Xij 



(fig f jo-} 



Xij ■ 1 



-i(cti —ctj ) 



(3) 



Atf = (/it/ji - klfn) - Ay • e^+^ . (4) 

The phase of the field \ is the compact U (1) gauge field, 
8ij = argXy — * 6ij + (da) tj with (da) y = aij - on, and 
the radial part, \ — \{Xij)\> is the so-called "resonat- 
ing valence bond" (RVB) coupling. The doubly-charged 
spinon-pair field A is an analog of a Cooper pair. 

At high temperature the RVB coupling is vanishing, 
X = 0, and the system is in the Mott insulator (or "poor 
metallic") phase. With decreasing temperature x be- 
comes non-zero, eventually enabling the formation of a 
spinon-pair condensate A = |(Ay)| and/or of a holon 
condensate b = {bi} @. Therefore, four phases |il Ifl Ho| 
may emerge: the Fermi liquid (FL) phase with 6^0, 
A = 0, the spin gap (SG) phase with b = 0, A ^ 0, the 
superconductor (SC) phase with b ^ 0, A ^ 0, and the 
strange metallic (SM) phase with 6 = 0, A = 0. 
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A compact Abelian two-Higgs model (cA2HM) in three 
dimensions with a U(l) gauge link field 8, a singly- 
charged holon field <£>i, and a doubly-charged spinon-pair 
field $2 is a feasible phenomenological model to describe 
the phase structure of the iwo-dimensional spin-charge 
separated quantum system 0,0. This model suffers from 
several limitations [2j to describe high-T c materials in the 
underdoped regime, since the model includes only Gaus- 
sian fluctuations about the mean field theory of the t-J 
model and does not reflect the SU(2) particle-hole sym- 
metry at half-filling. However, it seems plausible that 
the cA2HM can be applied in the less-studied overdoped 
regime (called sometimes the SM region 0, EI) where 
the particle-hole symmetry is explicitly broken. 

Since high-T c materials are type-2 superconductors, we 
restrict ourselves to the London limit in which the radial 
parts of both Higgs fields <&q — |<£>Q|e IVQ are frozen, 
|$i 2 1 = const. The action of the model is defined as 



Sa 



2HM 



-0/] cos Op 



Q=l I 



3(d<p Q + Q6) h (5) 



where Op is the standard lattice plaquette. The model (J5J 
obeys a lattice version of the £7(1) internal gauge sym- 
metry J3J|: 0^0 + da, <pq — > <pq — Qa. It describes the 
hole (electron) "constituents" by the dynamical holon <p\ 
and spinon if2 phases, which strongly interact via the dy- 
namical gauge field 0. The inverse gauge coupling (3 in 
Eq. JSJ, (3 = \f + Xb, is given by the diamagnetic suscep- 
tibilities Xf °f the spinon and Xb of the holon fields. The 
hopping parameters kq are connected to the doping con- 
centration x and the couplings t and J of as follows: 
K\ oc t ■ x and K2 oc J 

Multi-Higgs models similar to Eq. JSJ appear not only 
for theories of high-T c superconductivity but also in mod- 
els for metallic hydrogen [13 (using non-compact gauge 
fields an d Qa = —Q2 = 2) and for multi-band supercon- 
ductors 

The model (JSJ contains three kinds of topological de- 
fects: a monopole and two types of vortices, referred to 
as the holon and the spinon vortex d Q • The monopole 
has magnetic charge 27rfi, while the holon (spinon) vortex 
carries magnetic flux quanta 2irh (irk) of the gauge field 
0. Since the magnetic flux is conserved, one monopole 
is simultaneously a source of one holon vortex and two 
spinon vortices. The defects can be best studied in the 
dilute state characterizing the superconducting phase of 
the model with nonvanishing condensates. In this phase, 
the core of the holon (spinon) vortex is in the FL (SG) 
phase. The monopole core is in the SM phase where both 
condensates are zero. 

The structure of the three-dimensional phase diagram 
of the model is quite complicated. However, the faces 
and edges of the 3D phase "cube" (i.e. the limiting cases 
of vanishing or large couplings (3 and ^1,2), respectively, 
can be related to various well-known and sometimes non- 



trivial systems. These limiting cases are shown schemat- 
ically on the faces of the 3D phase cube presented in 
Fig. ^a) and are discussed in detail below. The interior 
of the cube in Fig.^a) is our qualitative prediction of the 
phase structure of the cA2HM. The shaded («i, K2) plane 
should be compared with Fig. ^b) which represents the 
result of our numerical investigation for (3 — 1.0. The nu- 
merical result is consistent with the prediction and will 
be presented later. 
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FIG. 1: (a) The qualitative 3D phase diagram and (b) its 
numerically obtained 2D cross-section at (3 — 1. 



On the K2 = face the compact Abelian Q = 1 Higgs 
model (cAHMq = i) emerges. There the holon condensate 
b and the holon vortices are controlled by the compact 
U(l) gauge field. The phase of the spinon-pair field is 
disordered which implies condensation of the spinon vor- 
tices and, therefore, vanishing A. In the (/3,Ki) plane 
two phases 

El El El exist: (i) the broken/Higgs 
phase with b ^ at large /3/large K\ corresponding to 
the FL phase, and (ii) the confining/symmetric phase 
with b w at small /3/small k\. The confining phase 
is the SM phase with some properties of the FL phase. 
The holon condensate is nowhere exactly vanishing be- 
cause the phases are connected analytically. The holon 
vortices are almost massless and therefore dense in the 
SM phase, and they become heavy and dilute in the FL 
phase. At the k,\ = edge the model reduces to the 
compact U(l) theory which is confining at any (3 due 
to the monopole plasma The plasma destroys the 

holon condensate, b = 0. At the edge (3 — > 00 the model 
becomes the 3D XY model with a second order transi- 
tion at ki.c — k^ y w 0.45420 . . . Since away from 
this limit a relevant global symmetry is absent, no sin- 
gularities in the thermodynamic quantities are expected. 
With decreasing (3, the SM and FL phases, being ana- 
lytically connected, can be thought to be separated by a 
Kertesz line [2(J, across which the holon condensate is a 
smooth quantity. 

On the face K2 — > 00 the model © reduces to a 3D 
XY-Z2 model with the action 



SxY- 



^2 P ap ~ ni ^2, m cos ( d( p)i 
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This model describes a XY-type matter field (p — ipi — 
(P2/2 which interacts via exchange of a 2^2 gauge field 
a 1 = ±1. The model - studied numerically in Ref. [2l| 

- contains three phases which, according to the adopted 
classification, are: (i) the SM phase at small /3/small K\, 
(ii) the SG phase at large /3/small K\ and (iii) the SC 
phase at large /3/large K\. 

For not too large K\ the SM and SG phases are sep- 
arated by a second order phase transition of the Ising 
universality class across a line (3 c {k\) with the 3D gauge 
Ising (gl) model limit H^: /3 C (0) = f3f w 0.7613 ... At 
large (3 the SG phase and SC phases are separated by 
a second order phase transition in the 3D XY univer- 
sality class, lim^oo /Ci >c (/3) = n^ Y . Inside the phase 
diagram these transition lines merge and continue as 
a single second order transition in the XY universal- 
ity class along a line ki iC (/3) which ends at the (3 — 
edge in ki jC (0) = k* y w 1.6. Here the model is re- 
duced to a (modified) XY model with the action Sx y = 

— J2i 1°S cosh[fti cos(dip)i]. 

The model © has exotic phases predicted for strongly 
correlated electron systems 23] and it is characterized 
by so-called "visons" - fractionally charged excitations 
related to the presence of Z2 degrees of freedom. A vi- 
son coincides with the spinon vortex at K2 — > 00, while 
the holon vortex turns into a XY vortex. The XY phase 
angle ip becomes the phase of the so-called "chargon" 
particle. The SG phase corresponds to the fractionalized 
phase where visons are absent and chargons are free par- 
ticles. In the SM phase the visons are condensed, and 
chargons are confined. The SC phase corresponds to a 
superfluid state where both visons and XY vortices are 
dilute and chargons arc free. 

On the «i = face the holon condensate vanishes, 
6 = 0, and the spinon-pair condensate A is described by 
the cAHMq = 2 model. The phase diagram in the (/3,K2) 
plane can be deduced from results of Ref. [II HI: The 
SM phase is located at small /3/small K2, while the SG 
phase is residing in the large /3/large K2 corner. The 
phases are separated by a second order phase transition 
line which ends at f3 c = /3f : in the large K2 limit and at 
«2,c = K^- Y in the large (3 limit. 

On the Ki — » 00 face the model (J5J) reduces to a XY 
model, Sxy = -^^i^i^^i, f = 2 <^i +^2, which 
controls the spinon condensate. Due to the constraint 
9 = Aipi + 2ttI, I G Z5, the holon vortices are suppressed 
and therefore b ^ in the whole {(3 ^2) plane. The phase 
diagram is divided by a straight line of a second order 
AY"-like transition parallel to the f3 axis at k,2,c{(3) = 
k* y ■ This line separates the FL phase (with condensed 
spinon vortices and A = 0) at K2 < from the SC 

phase (with suppressed spinon vortices and A / 0). 

On the (3 = face we obtain a two-Higgs system in- 
teracting ultra-locally via a non-propagating gauge field. 
After integrating out the gauge field in the Villain rep- 
resentation |2^], one gets a XY model with the effective 



AY-coupling, k^j- = kY K 2 / ( K Y + ^ K Y) where the su- 
perscript V indicates the Villain coupling. The phase 
diagram in the (ki,K2) plane contains two phases: the 
SC phase with non-zero condensates A and b in the large 
K\ /large K2 corner and a SM-FL phase in the remaining 
part of the phase diagram. The phases are separated by 
a second order AV-type phase transition which starts at 
Wi,c = K * Y at the K2 — » 00 edge and ends at K2, c — k% y 
for Ki — > 00. At these two edges the model is reduced to 
a modified and a usual XY model, respectively. The SM- 
FL phase is actually the FL phase at large ki, where the 
holon condensate is non-vanishing, and the SM phase 
at large K2, where the holon condensate tends to zero. 
Those phases are analytically connected, similar to the 
«2 = face considered above. 

On the (3 — > 00 face the system reduces to two de- 
coupled XY models describing the holon and spinon su- 
perfluidities. The phase diagram in (ki,K2) includes all 
SM, SG, FL, SC phases as shown on the upper face 
of Fig. nja). The predicted 3D phase diagram implies, 
in particular, that at sufficiently strong gauge coupling, 
[3 < (3f-, the SG phase ceases to exist. 

Performing Monte Carlo simulations we have investi- 
gated the phase diagram of the cA2HM. We made an 
exploratory first study on a 16 3 lattice, choosing three 
values of the gauge couplings, [3 = 1.0, 1.5 and 2.0, on a 
dense grid of points spanning the (k\,K2) hopping param- 
eter plane over the range < ki_2 ^ 2.0, 2.5. In all figures 
we show the results for the strong coupling (3 — 1.0 only. 

The density of the monopoles, p mon , is plotted in the 
upper panel of Fig. over the («i,/«2) plane. With in- 
creasing hopping parameters n\ or K2 the monopole den- 
sity gets suppressed. As shown in the middle panel of 

Fig. the density of the holon vortices p^ lt (spinon vor- 

(2) 

tices Pvo r t) significantly drops down with increasing K\ 

The connectivity C^ 2 -* = C^ 1 ' 2 ^ (i? max ) of the vortex 
clusters is derived from the cluster correlation functions 

tfinnjto E^S. QM&MS(R -x + v)) 

(0 c (x) = 1 or is the characteristic function of the vor- 
tex cluster c, i? m ax the maximal distance on the periodic 
lattice) gives a clear view of the phase diagram of the 
model [l7|, I20L |26J . Any Higgs condensate suppresses the 
proliferation of the corresponding vortices: they are pre- 
vented to percolate over infinitely long distances. We 
show the connectivities C IR ' for the spinon (holon) vor- 
tex clusters in the lower panel of Fig- El T ne correspond- 
ing phase diagram is shown in Fig. ^b) . It agrees with 
the cross-section through the 3D phase cube marked by 
the shaded plane in Fig.^a). We plotted letters A, B, C 
and D in Figs. H^a) and (b) to identify the same points. 

To clarify the nature of the phase transitions, we stud- 
ied the volume dependence of the average plaquette and 
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holon vortex spinon vortex 



FIG. 2: The monopole density (top), the vortex densities 
(middle) , and the vortex connectivities (bottom). Left: holon 
vortices; right: spinon vortices. 

both link contributions to the action JSJ as well as their 
susceptibilities, respectively, in different regions of the 
phase diagram at (3 = 1.0. Fig. (left) clearly exempli- 
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k, 2-state signals 

Holon link expectation value Histograms of action parts 
FIG. 3: Signatures for the first order transition. 

fies the first order nature of the transition (observed in 
all parts of action) by the jump developing in the holon 
link vs. K\ at fixed Ki in the crossing region of the tran- 
sition lines. Already at lattice size 32 3 we do not observe 
tunnelings for the selected K\ values within 5 x 10 5 it- 
erations. In Fig. |3 (right) we present typical two-state 
signals (for a 12 3 lattice) of all three terms in the action 
at («i, K2) near the transition. The two-state signal with 
respect to the plaquettc and the holon link becomes very 
weak when one goes to smaller K\ along the (except for 
the direct transition between the SM and the SC phases) 
horizontal dark-dotted (blue) line. 

It is known that for K\ —> the transition vs. K2 is 
of second order. Similarly, for large K\ the transition is 
most likely of second order, again in the XY universality 
class. We found that for the largest volumes 40 3 — 48 3 



under study the increase of the spinon link susceptibility 
at ki = 2.0 stops as function of the lattice size as ex- 
pected for the XY model at m — > 00 |2jJ. Selecting two 
«2 values outside the crossing regime, where one of the 
transition lines (the light-dotted [red] one) runs vertical, 
we observe that there is no thermodynamic transition vs. 
K\ for the smaller k 2 , in agreement with what could be 
anticipated from the limit K2 0. 

Summarizing, we predicted the qualitative 3-D phase 
diagram of the cA2HM proposed as an effective model 
for high-T c cuprates in the overdoped regime. At strong 
gauge coupling (however (3 > /Jf 1 ) we found two transi- 
tions associated with the pattern of vortex percolation 
and identified the Fermi liquid, spin gap, superconduc- 
tor and strange metallic phases in agreement with the 
proposed phase diagram. The percolation transitions are 
accompanied with thermodynamic phase transitions ex- 
cept for small holon hopping parameter. First hints for a 
changing order along the thermodynamic transition lines 
are found. In the region where the two transition lines 
merge, a direct phase transition between the SC and SM 
phases exists and is found to be first order. 
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